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Abstract 


The  problem  of  a  point  electric  dipole  moving  over  a  dispersive 
dielectric  half-space  is  studied.  The  dipole  is  located  in  the  free 
space  above  the  dielectric  and  is  assumed  to  be  time  harmonic  in  its 
rest  frame,  oriented  perpendicular  to  the  interface  and  moving  parallel 
to  it.  Previous  work  in  this  area  has  been  mostly  confined  to  the  use 
of  a  plane  wave  as  a  source  or  a  non-dispersive  dielectric. 

Solutions  for  this  problem  are  obtained  using  integral  transform 
techniques.  The  integration  is  performed  for  the  free  space  region  by 
making  an  asymptotic  expansion  for  one  integration  and  then  usintt  the 
saddle  point  approximation  for  the  remaining  integration.  The  solution 
obtained  in  this  manner  for  the  vector  potential  is  then  used  to  generate 
the  electric  and  magnetic  fields  in  both  the  rest  frame  of  the  source  and 
the  rest  frame  of  the  dielectric.  This  yields  the  reflected  radiation 
fields  for  the  case  of  an  arbitrary,  dispersive  dielectric  and  the  lateral 
wave  fields  and  surface  wave  fields  for  the  case  of  a  lossless  plasma. 

The  field  patterns  are  distorted  by  the  relative  motion  of  the  source 
and  dielectric.  In  the  rest  frame  of  the  source,  all  three  waves  exhibit 
the  frequency  of  the  source  and  in  the  rest  frame  of  the  dielectric,  they 
all  exhibit  some  form  of  doppler  shift  (different  from  the  primary  wave). 
The  criteria  for  existence  of  the  lateral  wave  is  not  modified  by  the 
relative  motion  of  the  source  and  dielectric;  but  the  criteria  for 
existence  of  the  surface  waves  is  greatly  modified  by  the  relative  motion. 
For  velocities  greater  than  some  critical  velocity,  a  new  type  of  surface 
wave  comes  into  existence. 
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1.  INTRODUCTION 


A  Moving,  temporally  dispersive  material  is  unusual  in  that  it  is 
both  temporally  and  specially  dispersive  as  observed  in  any  inertial 
reference  frame  other  than  its  rest  frame.  In  order  to  investigate  the 
interaction  of  electromagnetic  sources  with  such  a  material,  the  problem 
of  a  point  electric  dipole  moving  over  a  dispersive  dielectric  half-space 
has  been  studied.  The  dipole  is  located  in  the  free  space  above  the 
dielectric  end  is  assumed  to  be  time  harmonic  in  its  rest  frame,  oriented 
perpendicular  to  the  interface  and  moving  parallel  to  it. 

Previous  work  in  the  area  of  scattering  from  a  moving  half-space  has 
beer,  mostly  confined  to  the  use  of  aplane  wave  as  a  source  or  a 
non-dispersive  dielectric.  Particular  problems  involving  the  reflection 
of  plane  waves  from  a  moving  half-space  have  been  investigated  in  detail 
by  Tai  (1965),  Yeh  (1965,  1966,  1967)  Pyati  (1966)  and  others;  and  just 
recently  Shiozawa  and  Haxama  (1968)  have  solved  the  general  problem  of 
reflection  and  transmission  of  a  plane  electromagnetic  wave  at  the  inter¬ 
face  between  a  stationary  dielectric  medium  and  a  dielectric  medium  moving 
in  an  arbitrary  direction  parallel  to  the  interface.  It  appears  that  the 
only  work  that  has  been  done  on  the  scattering  of  dipole  radiation  by  a 
moving  half-space  was  done  by  Pyati  in  1966.  He  investigated  the  problem 
of  a  vertical  or  horizontal  dipole  positioned  over  a  moving,  dielectric 
(non-dispersive)  half-space  and  solved  for  the  reflected  radiation  fields 
in  the  rest  frame  of  the  source.  No  effort  was  made  to  investigate  the 
effect  of  the  motion  on  the  lateral  wave  or  surface  waves  that  might  be 
present. 

In  the  work  that  follows,  we  utilise  two  inertial  reference  frames 
designated  S  and  S’  where  it  is  assumed  that  reference  frame  S'  is  moving 
with  velocity  v  in  the  i  direction  with  respect  to  reference  frame  S  and 
that  they  coincide  at  t  =  0.  The  dielectric  half-space  is  stationary  in 
reference  frame  S  and  occupies  the  region  x  <  0.  The  dipole  is  stationary 
in  reference  frame  S'  and  is  positioned  at  the  point  (H,  0,  0)  in  that 
frame.  The  region  x  >  0,  in  which  the  dipole  is  located,  is  assumed  to  be 
free  space. 

Solutions  for  this  problem  are  obtained  using  integral  transform 
techniques.  Section  2  develops  integral  formulations,  involving  three 
integrations,  for  the  vector  potential  in  reference  frame  S  for  both 
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regions  of  spec*  and  carries  out  one  of  the  integrations.  Zt  appears 
that  this  is  as  far  as  the  work  can  be  carried  out  in  general.  3eyond 
this  point  sqm  sort  of  approximation  aust  be  made.  This  is  done  in 
sections  3  through  6  where  solutions  are  found  for  the  fields  '-t 
observation  points  far  from  the  image  source. 

Zn  section  3  the  original  contour  for  the  integral  is  converted  to 
contours  around  each  singularity  in  the  integrand*  and  a  change  of  variable 
is  made  to  gain  control  over  the  asymptotic  expansion  variable  that  is 
subsequently  used.  These  singularities  consist  of  two  branch  points  and 
one  or  more  poles.  The  integrations  along  the  branch  cuts  are  carried  out 
in  sections  4  and  5  and  yield  the  reflected  radiation  fields  for  a 
generalised  dispersive  dielectric  and  the  lateral  wave  fields  for  a 
lossless  plasms  dielectric.  The  integration  around  the  poles  is  carried 
out  in  section  6  and  yields  surface  waves  for  the  case  of  a  lossless 
plasma  dielectric.  Zn  each  case,  the  saddle  point  approximation  is  used 
for  the  final  integration.  The  vector  potential  found  in  this  manner  is 
then  used  to  generate  the  electric  and  magnetic  fields  (E  and  I)  in 
reference  frme  S.  These  fields  are  then  transformed  to  reference  frame 
S'  to  obtain  the  electric  and  magnetic  fields  (S'  and  S')  in  the  rest  frame 
of  the  dipole. 

2.  INTEGRAL  FORMULATION 

Zn  this  stetion,  an  integral  formulation  for  the  vector  potential  in 
the  rest  frame  of  the  dielectric  is  developed.  Since  the  problem  that  is 
being  considered  is  homogeneous  with  respect  to  the  extension  space 
coordinates  y,  s  and  t,  it  lends  itself  nicely  to  the  use  of  Fourier 
transform  techniques  for  the  solution  with  respect  to  these  coordinates. 
This  leaves  us  with  a  one -dimensional  boundary  value  problem  to  aolve 
similar  to  the  familiar  Scmmsrfsld  Half-Space  problem.  However,  for  the 
problem  we  are  considering  the  source  is  men's  interesting  sines  it  is 
both  temporally  and  specially  dispersive  due  to  the  notion  of  the  source. 

Before  we  can  solve  the  boundary  value  problem  in  reference  frame  S, 
we  must  first  obtain  a  formulation  for  the  source  as  it  would  he  observed 
in  that  reference  frame.  Zn  reference  frame  S',  the  source  is  described 
by 
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jw  t* 

J’(r',t*)  *  a  I#  5(x‘-H)4(y* )«(*')•  *  (1) 

We  would  like  to  convert  3* (?' ,t * )  to  fnr^ut^cy  space  so  that  we  nay 
conveniently  transform  it  to  the  rest  fraate  of  the  dielectric.  The 
four-dimensional  Fourier  transform  of  (1?  X/;  given  by 

5MK 

Si  I  e  A  ) 

3'(k’,«  )  * - ! — jjj — * - *-  (2) 

where  the  kernel 

-jtf’t*  +  jit**?' 

S? 

has  been  used  in  taking  the  transform.  We  will  now  use  the  relativistic 
spectral  transformations  (Holmes,  1968)  to  transform  the  spectral  form 
of  the  current  density  given  by  (2)  to  reference  frame  S.  The  transfor¬ 
mations  are  given  by 


Jx(k,W)  a  J^k’.tt*) 

(3) 

Jy(k,«)  = 

00 

J <£,«)  =  r[J’(k’,«»)  ♦  vp’(Ef,w’)3 

2>  Z 

(5) 

The  charge  density  P' (£’  ,«* )  can  be  obtained  from  3* (k*  »w’ )  by  using  the 
continuity  equation  converted  to  frequency  space.  The  spectral  form  for 
the  current  density  in  frame  S  then  is  givsn  by 

I  jk,H  «  vk, 

J(k,w)  *  e  6(m-(-yn  *  vkj))C*  ($) 

3 

where  the  transformations  for  the  components  of  the  wave  four-vector 
(Jackson,  1962)  have  been  used.  Now  inverting  (6)  with  respect  to 
we  obtain 
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^(XikgtkjtW) 


Z  W 

— *-  ♦  vk,)>Cft  «<x-H)  ♦ 

r  ft*  ~  3 


(7) 


which  is  ths  dssirsd  source  tens. 

The  one -dimensional  boundary  value  problem  we  must  solve  is  illustrated 
in  figure  1  where  use  has  been  made  of  (7)  and  the  fact  that  in  free  sj>ace 
the  fore  of  the  differential  equation  for  the  vector  potential  is  invariant. 
Since  the  solution  of  this  type  of  problea  in  the  manner  we  are  proposing 
is  well  known,  the  details  have  not  been  included1.  However,  we  would 
like  to  make  a  few  specific  comments.  At  first  glance  it  appears  that 
difficulty  sight  be  encountered  in  finding  the  particular  solution  to  the 
differential  equation  in  region  1.  However,  due  to  the  delta  functions  in 
the  source  terms  it  can  be  readily  generated  using  the  free  space  Green's 
function  for  the  differential  equation.  Also,  it  turns  out  that  the 
vector  potential  obtained  from  the  solution  of  this  boundary  value 
problea  contains  the  factor 


vk3)) 


so  that  the  inversion  with  respect  to  w  can  be  performed  exactly.  The 
vector  potentials  correponding  to  the  scattered  and  transmitted  field 
are 


AlzS(r’t) 


jvrp  i  3r“,(t  -  *> 

* - S-i-e  c 

2w#(2irr 


•  K-h  -jf 

dk33dk2 

•w  m  m 


(8) 


-JU  I 

S*c*r»t)  *  I 

1X5  2<2w)2 


jr«#(t  -  %  Z)  . 


fit. 


2  K(k33  ♦ 

2rv  ■  33  c2  ~r 


(KCr+h)(X+h) 


-jfi 


1.  Considerable  detail  is  included  in  the  original  work  (Holmes,  1968) 
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A2zT  * 


j u  i  v  jrw,(t  ‘  *T  z)  -  K  -jf 

e  •  _  c  //  J 


u  (210' 
o 


TK+hT 


e  dk.  dk_ 
43  2 


(10) 


and 


■ill,  3r“."  '  -  r2  5‘-cr 

^  _ •  •  -  c  ff  fl  V 


2xT 


(210' 


II  ^ 


VUl 

-l)(k33  ♦  -jt) 
c 


U  (  Ke+h  )  {  K+h ) 


r  -  -  - 


♦  rrr^r-]  e  ^  2  dk^dk. 


K€  +h 
— r  - 


‘33  2 


(11) 


where 


cj 


K  -  _  k2  k2  vi/2 

-  (c2  k2  k33} 


(12) 


„  r  9  ?  ^  9  o  wo 

!>  =  tfr  T  ♦  ',k33)*  -  r  <k33  *  y->  -  k2] 


(13) 


0  *  c 

-r  rjw  =  T(u>  ♦  vk,_) 

$  33 


(14) 


fx  *  K(K+x)  +  k33  T(z-vt)  +  k  y 


(15) 


^2  *  ^  **  hx  +  k33  F(z-vt)  +  k2y 


and  the  change  of  variable 


k 


3 


r(k 


33 


vw 


(16) 
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has  been  made  to  simplify  the  expressions.  It  appeal's  that  this  is  as 
far  as  the  work  can  be  carried  out  in  general  and  that  some  sort  of 
approximation  must  be  made  in  order  to  perform  the  remaining  two 
integrations.  This  is  done  for  region  1  in  the  following  sections. 

3.  SCATTERED  FIELDS 

He  would  like  to  evaluate  the  integrals  for  the  scattered  fields  by 
obtaining  the  first  term  of  an  asymptotic  expansion  for  the  integration 
with  respect  to  k33  and  then  using  the  saddle  point  approximation  for  the 
integration  with  respect  to  We  can  gain  control  of  the  asymptotic 
expansion  variable  by  putting  the  Integral  into  the  form  (Van  Der  Waerden, 
1950) 

I  =  /  P(u)e“Xu  du  (17) 

A  convenient  choice  for  X  is 

w 

X  *  r'  -2-  (18) 

-  c 


where 

r'  =  [(H+x)2  t  2(z-vt)2  +  y2]1/2  (19) 

and  is  the  distance  in  wave  lengths  from  the  image  source  to  the  field 
point  in  reference  frame  S’.  Then  for  some  suitably  large  X,  the  first 
term  in  the  asymptotic  expansion  should  in  most  cases  represent  the 
integral  adequately. 

In  the  u  plane,  the  integrands  under  consideration  have  singularities 
in  the  form  of  branch  points  and  poles.  The  original  contour  can  of  course 
by  the  Cauchy-Goursat  Theorem  be  replaced  by  paths  around  these  singu¬ 
larities  plus  a  suitable  closure  path.  The  closure  path  will  be  chosen  as 
the  right-half  plane,  to  ensure  that  there  will  be  no  contribution  to  the 
integral  from  the  closure.  The  branch  cuts  will  be  taken  from  the  branch 
point  along  a  path  parallel  to  the  real  axis.  This  forces  the  integrand  to 
decrease  nnst  rapidly  along  the  branch  cut  paths  as  you  move  away  from  the 
branch  point  and  will  aid  in  our  evaluation  of  the  branch  cut  integrations. 
Our  original  integral  then  becomes 
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I  =  -/  P(u)e"Xu  du+2*j  l  R.  SIGH(z-vt)  (20) 

i  1 

Branch  cut  paths 

where  is  the  residue  at  the  ith  pole  and  SIGN(z-vt)  is  11  depending  on 
the  sign  of  the  quantity  (z-vt). 

Let  us  now  consider  the  procedure  that  will  be  followed  in  evaluating 
the  branch  cut  integrations.  The  first  step  will  be  to  transform  the 
branch  point  to  the  origin  and  then  to  utilize  the  second  branch,  P(u~), 
fcr  the  path  below  the  cut  and  the  first  branch,  P(u+),  for  the  path  above 
the  cut.  The  integration  is  completed,  by  expanding  the  integrand  in  a 
series  and  then  keeping  only  the  first  term.  The  resultant  expression  for 
each  branch  cut  integration  becomes 


QO 

U  =  -SIGN(z-vt)e~Xd  B(k9)  /  nN  e-Xn 


where  N  and  B(k2>  depend  on  the  integrand  under  consideration  and  d  is  the 
location  of  the  branch  point.  Evaluating  the  integral,  this  becomes 


IB  *  -SIGN(z-vt)e~Xd  B(k2)A’(N+1)r(H+l) 


The  integration  with  respect  to  k^  will  be  done  by  a  saddle  point 
approximation.  For  the  branch  cut  contributions  the  integral  will  be 


of  the  form 


Ag  *  -SIGN(z-w)  /  BB(k2)e  *  dkj 

A  '  '  “°° 


f(k2) 


where 


f(k2>  =  jCTw  (t - 2  z)  -  k2y]  -  Ad 

c 


The  saddle  point  approximation  gives  us 


ihunWibUnnuiiW' > ~ ... ..u«. W  iji  uttphit dd 1  Cm*' 
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AB  = 


-j  SIGN(z-vt)  /Tn 


r(K+l)BB(k20)e 


-f<k20> 


,<N.1>/S  f(k20) 
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3k! 


(25) 


where  the  saddle  point  k2Q  is  found  from 


(26) 


The  integration  with  respect  to  k for  the  pole  contributions  to  the 
integral  will  also  be  done  by  the  saddle  point  approximation.  The 
details  are  discussed  in  section  6. 

The  above  work  describes  in  general  the  technique  that  will  be  used 
to  evaluate  the  integrals.  The  actual  integrals  will  be  evaluated  in 
sections  4  througn  6  using  the  above  results.  We  would  now  like  to 
investigate  the  nature  and  location  of  the  singularities  that  we  will 
have  to  deal  with  in  those  sections. 

Before  examining  the  integrands  for  singularities  we  must  transform 
them  to  the  u  plane.  They  will  be  given  by 


P(u) 


W(u) 


(27) 


where  W(u)  is  the  integrand  in  the  plane  expressed  as  a  function  of  u. 

^33 

We  therefore  need  to  find  k^,  K»  and  h  as  a  function  of  u.  Using 

(9),  (15)  and  (17)  we  find  that  the  first  three  are  given  by 


b> 


33 


=  -r ' 


2  2 

(rf  -y  ; 


CjuHz-vt)  +  (Htx)SIGN(z-vt) 


k2c2  2 
u2+(l-  -i—Ml-  JL^)] 

u>  r» 

o 


(28) 
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and 


r'« 


c(r 


V-**- . c-ju(H+x)  ♦  Jr(*-vt) 

»  -y  ) 


- 2T - 

KC  2 

U2  +  (1.  X-Xl-  Xlr)]  (30) 
eT  r' 


The  function  h(u)  cannot  be  found  unless  the  dielectric  has  been  specified. 
For  the  case  of  a  lossless  plasea 

w2 

t  *  1  -  (31) 

r  w2 

and  h  takes  on  the  particularly  sieple  for* 


h  a 


(32) 


We  ore  now  ready  to  exaaine  the  integrands  for  the  pertinent  singular 
points.  Each  singular  point  of  P(u)  that  is  enclosed  by  the  controur  will 
make  a  contribution  to  the  integral.  In  lieu  of  directly  finding  which 
points  are  enclosed  by  the  controur,  and  then  evaluating  their  contribution 
to  the  integral,  we  have  exanined  the  contribution  associated  with  each 
singular  point  and  then  retained  only  those  contributions  :*hich  correspond 
to  waves  of  the  outgoing  type.  The  two  approaches  of  course  are  equivalent; 
however,  the  latter  technique  is  more  tenable  when  a  aultiple  integral  is 
involved,  since  the  location  of  the  singular  points  will  generally  £•«  a 
function  of  the  subsequent  integration  variable. 
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From  examining  (8)  and  (9)  we  see  that  the  integrands  under  consideration 
have  branch  points  where 

K  =  0  (33) 

and 


h  =  0 


(34) 


and  that  the  integrand  in  (9)  may  have  one  or  more  poles  where 


K  e  +  h  =  0 
-  -r 


(35) 


In  the  u  plane,  the  branch  points  corresponding  to  (33)  and  (34)  are 
located  at 


2  2 

k,e*  2 

d  =  j  /  (1  -  -Ar-Hl  -  2-5-) 

u)  r* 

# 


(36) 


and 


d  *  ♦  I >■<*-«>  I 


r  ~2  2 
rU)  -u 


— £  -  k2] 
2  2J 


(37) 


The  poles  associated  with  equation  (35)  will  be  discussed  in  section  6. 

Our  ultimate  goal  is  to  find  the  scattered  electric  and  magnetic  fields 
in  reference  frame  S  (E  and  B)  and  reference  frame  S'  (E*  and  §*).  Due  to 
the  nature  of  the  integration  techniques  used  in  the  above  work  and  in  the 
work  that  follows,  these  fields  can  be  found  from  the  vector  potential  in  a 
simple  manner.  In  x,  k^»  k^,  u  space,  the  vector  potential  in  reference 
frame  S  has  the  form 

A  =  A(k2,k3,w)e  s  A(x,k2»k3,w)  (38) 

In  order  to  find  the  electric  fields,  we  need  to  find  an  expression  for  the 
scalar  potential.  The  scalar  potential  can  be  found  using  the  Lorentz  Gage 
condition  (which  is  invariant  in  free  space)  transformed  to  the  same 
coordinate  space  as  A.  The  scalar  potential  then  is  given  by 

c2 

♦(x,k2,k3,w)  *  ~  [KAx  t  k3Az3  (39) 


n 


f 


Using  the  vector  potential  and  (39),  the  electric  and  magnetic  fields  can 
be  generated.  This  -form  for  the  fields  would  in  general  have  to  be 
inverted  with  respect  to  k2,  and  w.  However,  in  the  above  work  and 
in  the  work  that  follows,  the  integration  process  is  such  that  the 
amplitude  portions  of  the  integrand  are  evaluated  at  some  particular 
w,  k2  and  k3«  It  is  therefore  not  necessary  to  go  through  this  integration 
process  again.  We  can  use  the  expressions  already  obtained  for  A_  and  A 
and  replace  w,  k^,  and  k2  in  the  expressions  for  the  fields  by  the  values 
used  in  the  particular  integration  process.  Proceeding  in  this  manner, 
the  fields  in  reference  frame  S  are  given  by 


where  k2  and  k33  depend  on  the  particular  singularity  under  consideration 
and  are  yet  to  be  determined. 
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4.  REFLECTED  RADIATION  FIELDS  FOR  A  LOSSLESS  PLASKA  DIELECTRIC 

The  generalized  technique  formulated  in  section  3  for  evaluating  the 
branch  point  contributions  to  Aix#  and  A^Js  will  be  applied  in  this  section 
to  the  branch  point  given  by  (36).  It  will  be  found  that  the  contribution 
fro«  this  branch  point  corresponds  to  the  reflected  radiation  fields. 

In  ordsr  to  use  (25)  which  gives  a  generalized  solution  for  the  branch 
cut  integration  we  need  to  find  the  saddle  point  and  to  calculate 

— j  (k2Q),  f(k2Q)*  N  and  BB(k2Q).  Performing  these  calculations  and  using 
3k2 

the  results  in  (25)  we  find  that  the  contributions  to  the  vector  potential 
fro*  the  branch  cut  integration  corresponding  to  equation  (36)  are  given  by 


hzS{*'x)  * 


u  T  K  K  K  W*  '  z>  -  *~3 

M  1  „  V  (H+X)  *  C2 

*  r  — Jr*  Ly-Vv-  3  e 


4ir 


(46) 


and 


*axS<?.*>  *  t2  r2 1 K 


(Hz-Vt)  V.  ^ 

1  F  cn-r# 


,2  12  c  -«  7k  V  +h  TTk  +h  ) 

r  -*-r»  *»  -i  -i 


lit  -h  )  1".cr(t ' 7  z)  - 

-*-r*  *»  _  c 

*  rTTcT+hr  • 


(47) 


where 


K  =  (Hfx) 
”  c  r' 


(48) 


h  *  4  Ce  r2A2  -  r2D2  -  *lr31/2 


-•  c  -re 


(49) 


-re 


*  6_ 


iw  * 


Tw  A 

o 


(50) 
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A  *  1  ♦ 


3r(z-vt) 


D  s  rS.2;vfl  t  B  (52) 

r 

Equations  (46)  and  (47)  correspond  to  the  reflected  radiation  fields. 

We  observe  that  the  fields  go  down  as  1/r*  where  r'  is  the  distance  from 
the  image  source  to  the  field  point  and  that  although  the  frequency 
distribution  is  the  same  as  that  of  the  primary  source*  the  distribution 
is  centered  on  the  image  source.  This  means  that  whenever  the  velocity 
is  non-zero,  the  frequency  of  the  reflected  radiation  field  as  observed 
in  reference  frame  S  will  differ  from  that  of  the  primary  field  everywhere 
except  on  the  interface  and  a  vertical  axis  through  the  dipole.  The 
frequency  distribition  is  given  by 

~  *  r[l  +  ~  (53) 

iti  c  r* 

D 

We  also  observe  that  in  the  rest  frame  of  the  source,  the  reflected  radiation 
field  will  not  have  a  doppler  shift.  It  should  be  noted  that  (46)  and  (47) 
are  valid  for  an  arbitrary  dispersive  dielectric.  When  v  *  0,  (46)  and  (47) 
reduce  to  the  well  known  result  for  the  zero  velocity  case  (Banos,  1966, 
pg.  35  6  179). 

{text  we  would  like  to  find  the  electric  and  magnetic  fields  in  reference 
frame  S  and  reference  frame  S’ .  Equations  (40)  through  (45)  give  these 
fields  in  reference  frame  S  in  terms  of  the  values  of  and  k 2  used  in  the 
above  integration  process  and  given  here. 


k  =  -A  X 

2  c  F 


w  -  •  ru-vt) 
33  *  c  F 


Using  these  In  (40)  through  (45),  the  electric  and  magnetic  fields  in 
reference  frame  S  are  given  by 


I 


. . . . 
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"xR 


r»*2/,  V  ru-vt).2  (H+x)2,* 

■  -  V-  rc.'-Vt-y; [r  (1  ♦  ?-?— >  -7*-** 

6  P*  • 


(Ifc*)  ../Hx-vt)  .  vw  , 
-p-rt-p — ♦  c)axr3 


(56) 


P  _  ^W»y _  r(H+X)  k  r/r(z-Vt)  V. 
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-  c  r 1 


(57) 
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xR 


(59) 


■  rr^(l*yt)  .  Vi.  i 

AxR  **  r^~P  +  c3AxR3 


(60) 


and 


(d 


B*r‘ 


j  *f£AxR 


(61) 


where  AxR  and  A^  are  tha  x  and  z  components  of  tha  vector  potential  for 
the  reflected  radiation  field  given  by  (46)  and  (47).  He  would  now  like 
to  tranafcra  the  electric  and  aagnetic  fields  for  the  reflected  radiation 
field  to  reference  frame  S'.  The  transformation*  are  well  known  (Jackson* 
1966)  and  yield 


E» 

xR 


*  * 1  r» crw  ■ 


)A 


xR 


(H+x») 


**AzR3 


(62) 


r.  -  ,  V  rr(Htx')  *’  *  •, 

V  •  3  TFTa  - AxR  +  F  zRJ 
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^zR  '  3  Ta  t~F (p-  +  c)AxR  ■  (1  "  ^,A^] 


,2  zR 


Bt  -  ~jV  r  v  (H»x* )  +A  1 
xR  fcAr '  ~  w®  ’D 


c  r’  xR  zR 
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V  '  fST  L'r  {p-  A  *  c  ■  72"  "  )AxR  +  ~P - AzR3 


BzR  =  3  c  r^*  AxR 


where  A  s  X  ♦ - -r  . 

c  r 

The  reflected  radiation  fields  for  an  arbitrary*  dispersive  dielectric 
are  given  by  (5o)  through  (61)  in  the  rest  frame  of  the  dielectric  and  by 
(62)  through  (67)  in  the  rest  frame  of  the  source.  Although  these 
equations  are  not  terribly  complex,  it  is  nevertheless  difficult  to 
determine  the  effect  of  the  motion  on  the  fields  by  merely  inspecting  the 

2  -  -  iiflpA  ,  ,  41fr  « 

equations.  Therefore,  the  quantities  (E„ •£*)(-——=—)  and  (E' »E '*)(-— -~-y—) 

K  H  U)  y  X  K  K  tv  M  1 

2  #  o  •  •  o  o 

u> 

were  plotted  for  a  lossless  plasma  with  X  *  *  5  in  order  to  observe 

some  of  the  effects  caused  by  the  relative  motion  cf  the  source  and  dielectric. 
Figure  2  shows  the  field  patterns  in  the  x-z  plane  for  x  >  0  in  a  polar  plot. 
For  S  =  0.0,  the  pattern  is  a  symmetrical  two  lobe  pattern.  In  the  rest 
frame  of  the  dielectric,  as  the  velocity  is  increased,  the  lobe  behind  the 
source  shrinks  and  the  lobe  in  front  on  the  source  "peaks  up"  and  bends 
down  toward  the  direction  of  motion.  This  is  probably  what  should  be  expected 

2.  R*  is  an  arbitrary  fixed  distance  from  the  image  source  as  measured  in 
reference  frame  S. 
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since  the  reflected  wave  is  highly  dependent  on  the  primary  wave,  and 
the  primary  wave  has  this  same  behavior  (Fujioka,  1966).  In  the  rest 
frame  of  the  source,  as  the  velocity  of  the  dielectric  increases  the 
lobe  in  the  direction  of  motion  of  the  dielectric  grows  larger  and  bends 
down  toward  the  interface  and  the  other  lobe  shrinks. 


5.  LATERAL  WAVE  FIELDS  FOR  A  LOSSLESS  PLASMA  DIELECTRIC 


The  generalized  technique  formulated  in  section  3  for  evaluating  the 
branch  point  contributions  to  A^xs  and  A^s  will  be  applied  in  this  section 
to  the  branch  point  given  by  (37).  It  will  be  found  that  the  contribution 
from  this  branch  point  corresponds  to  the  lateral  wave  fields. 

Proceeding  as  in  section  4  we  find  that  the  contributions  to  the 
vector  potential  from  this  branch  point  are  given  by 


.  .  (1-X)1/2  >I/4  ef<IC20> 

lzL  ^  R' 2l/x  -  i—ip-  /T- jf/2 


(68) 


and 


A,  .  * 
lxL 


-jy  I  cU-X)1/2e  (  20>  T2  J  /hJ  t  -He2  -  1)  +  e 

s  > _ c  R  c  -re _ 

Irf 


2*0)  X 


1/4  R*2  e2 
-r# 


(H+x)  fiZxjS/T 


(69) 


where 


f(k20)  *  j«#cr(t  -  *)  -  /%  -  ~~  ]  (7o) 

C 


and 


5r* 


*  1 


r2(i 


+  — 
c 


V  TTz-VtT 


/TJ)4 


(71) 


The  above  is  for  the  case  when  the  dielectric  is  a  lossless  plasma.  We 
should  note  that,  since  the  location  of  the  branch  point  that  corresponds 
to  the  lateral  wave  is  independent  of  the  velocity,  the  criteria  for 
existence  of  the  lateral  wave  is  the  same  as  for  the  zero  velocity  case 
(Brekhovskikh,  1960).  When  v  =  0,  (68)  and  (69)  reduce  to  the  well  known 
result  for  the  zero  velocity  case  (Brekhovskikh  pg.  275,  Banos  pg,  184). 
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The  frequency  distribution  in  reference  frame  S  for  this  lateral  wave 
is  given  by 


!k,  u.idersspa* 


where 


S’  =  Cy2  +  r2(z-vt)2]1/2 


We  also  observe  that  in  the  rest  frame  of  the  source,  the  lateral  wave 
fields  will  not  have  a  Doppler  shift.  The  frequency  distributions  for 
the  lateral  wave  fields  (X  *  0.5)  and  the  reflected  wave  fields  are  shown 
in  a  polar  flot  in  figure  3.  The  plots  correspond  to  what  would  be  observed 
on  a  y-z  plane  intersecting  the  image  source.  It  should  be  noted  that  in 
making  the  lateral  wave  plot,  the  fact  that  the  group  velocity  for  the 
wave  is  not  the  speed  of  light,  was  properly  accounted  for  in  making  the 
ray  angle  transformations. 

We  would  now  like  to  find  the  electric  and  magnetic  fields  in  reference 
frame  S  (£  and  H)  and  reference  frame  S’  (E*  and  K' )  corresponding  to  the 
lateral  wave.  A  good  starting  point  is  (40)  through  (45)  which  give  these 
fields  in  terms  of  the  values  of  k33  and  k^  used  in  the  above  integration 
process  and  given  here. 

,,  .  JL  EiiSfll  /PJ  (7.) 


Using  (74)  and  (75)  in  (40)  through  (45)  we  obtain 

B  _  *  -  j  -1  X  /TJ  A.  r 

xL  J  c  R*  I2L 


V  *  1  T  *ui  '  r  r  ai«l] 


zL  J  c  R  lxL 
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F  i 
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and 


K 


bJl  s  rr cr  F  *■ ki«L  - (1  -  Z2  (1-X»W 


,2 


(89) 


The  lateral  wave  fields  for  a  lossless  plasma  are  given  by  (76) 
through  (81)  in  the  rest  frame  of  the  dielectric  and  by  (94)  through 
(89)  in  the  rest  frame  of  the  source.  In  analogy  with  the  discussion 
on  the  reflected  radiation  fields,  we  would  like  to  plot  quantities 
2  2 

like  (£'  *Ef )(~~ — )2  and  (E*  These  are  shown  in  figure 

l  l  y  1  c  L  y  1  c 

oo  oo 

4,  where  again  due  account  has  made  in  the  reference  frame  S  plots  for 

the  fact  that  the  group  velocity  for  the  lateral  wave  is  not  the  speed 
3 

of  light  .  The  above  quantities  are  plotted  in  a  polar  p)'  ->r  X  =  0.5, 

the  y-z  plane  and  (Htx)  -  0.  For  8  =  0.0,  the  patterns  l  •.  tropic. 

In  the  rest  frame  of  the  source,  as  the  velocity  is  increased.,  the  pattern 
becomes  elongated  in  the  direction  the  dielectric  is  moving.  It  almost 
looks  as  if  the  energy  being  radiated  by  the  source  is  being  dragged  along 
by  the  dielectric.  In  the  rest  frame  of  the  dielectric,  as  the  velocity 
increases  the  pattern  is  elongated  in  the  direction  of  motion  of  the  source. 

6.  SURFACE  WAVE  FIELDS  FOR  A  LOSSLESS  PLASMA  DIELECTRIC 

We  will  now  examine  the  contribution  to  the  integral  for  Aixs  for  a 
lossless  plasma  dielectric  due  to  the  pole  associated  with  the  function 
K  er  +  h.  In  order  to  find  the  location  of  the  pole,  we  must  solve  the 
equation 


u>2  2 

(1'ke  '  k2)1/2(1  “  -2-2-* - 2)  +  (i_k3  *  k2  _  "!)i/2  *  0  (90) 

"•  '2  rVu+sk.r  ~3  "2  w 

0-3  o 


where 


and 


k  ~  ■  —  v 
-3  u  33 
o 


k  *  —  k 
-2  u  k2 
o 


(91) 


(92) 


3.  See  Holmes  (1968)  for  a  more  complete  discussion  of  this  aspect  of  the 
problem. 
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for  the  value  of  k3  which  utiifiis  it.  Sine*  this  is  virtually  impossible 
to  do  analytically,  ws  will  assume  that  (90)  has  been  solved  and  continue 
with  the  evaluation  of  the  integral.  Then  later,  after  the  solution  has 
been  formulated,  (90)  can  be  solved  numerically. 

The  integration  with  respect  to  k33  is  then  just  2  jirR^  SIGN(x-vt) 
where  R^  is  the  residue  at  the  pole.  The  contribution  by  the  pole  to  the 
integral  for  Alx#  then  is  given  by 

U  I  .  Jr“.(t  -  T-*’  -  jr'  Iff, 

A1m.  ‘  -fee-1  [m  Rp  *  H  <"> 


where 


(94) 
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ip  *  <1  -  ijp  -  i2  *  ^>1/J  <97> 

)  (98) 

k^  equals  the  value  of  K3  which  satisfies  (90)  and  the  remaining  variable 
of  integration  has  been  changed  to  k^.  Using  the  saddle  point  technique 
to  perform  the  integration  with  respect  to  k2,  A^  becomes 
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where 


(99) 


po 


Pi 


k  =  k 
-2  -20 


and  k^0  is  obtained  from  the  solution  of 
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dk_  dK 

Equation  (101)  contains  the  functions  ^~-P  and  ^2-  .  These  can  be 

~2  -2 

obtained  from  (96)  and  (90).  Equation  (101),  from  which  the  saddle  point 
k2Q  is  found,  then  becomes 
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The  solution  for  Aixsg,  given  by  (99)  "goes  down  as"  (r*)  which 
implies  that  Aixsp  corresponds  to  a  surface  wave.  In  order  to  find  out  if 
a  surface  wave  is  possible,  the  phase  function  for  A^xjg  must  be  examined. 
When  v  *  0,  (99)  reduces  to  the  well  known  result  for  the  zero  velocity 
case  (Banos  pg.  124). 

A  parameter  associated  with  the  surface  wave  which  is  of  interest 
is  the  normalized  frequency  given  by 
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3  _  Tv  3 

3t  ■  p-  ot 


(104) 


He  note  that  in  the  rest  frame  of  the  source,  the  surface  wave  does  not 
exhibit  a  Doppler  shift. 

He  would  now  like  to  find  the  electric  and  magnetic  fields  in 
reference  frame  S  (2  and  H)  and  reference  frame  S'  (S'  and  H* )  corresponding 
to  the  surface  wave(s).  Proceeding  as  before  and  using  (91)  and  (92)  yields 
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The  surface  wave  fields  for  a  loss.. ess  plasma  are  given  by  (105) 

4 

through  (110)  for  the  rest  frame  of  the  dielectric  and  by  (111)  through 
(116)  in  the  rest  frame  of  the  source.  These  equations  are  not  explicit 
functions  of  the  coordinates,  but  are  expressed  as  functions  of  the 
coordinates,  but  are  expressed  as  functions  of  two  parameters  which  are 
obtained  from  the  simultaneous  solution  of  (90)  and  (102).  The  detailed 
computer  calculations  necessary  to  obtain  a  plot  of  the  surface  wave 
fields  are  in  progress.  However  a  preliminary  study  of  (90)  and  (102) 
has  been  performed  in  an  effort  to  find  how  the  surface  wave  modes  behave 
as  the  velocity  increases.  It  was  found  that  the  criteria  for  existence 
of  the  ordinary  surface  wave  is  modified  by  the  relative  motion  of  the 
source  and  dielectric.  As  the  relative  velocity  is  increased,  the  ordinary 
surface  wave  ceases  to  be  excited  in  all  directions  along  the  surface. 

The  wave  becomes  restricted  to  a  wedge  shaped  region  (0  5  wedge  angle  <2 ») 
behind  the  source.  This  is  illustrated  in  figure  5  by  the  top  three  curves 
which  show  the  wedge  angle  versus  velocity. 
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In  addition,  for  velocities  greater  than  some  critical  velocity,  a 
new  type  wave  comes  into  existence.  This  wave  is  a  surface  type  wave 
and  radiates  into  a  wedge  shaped  region  behind  the  source.  As  the 
velocity  increases,  the  wedge  angle  increases  from  zero  degrees  at  the 
critical  velocity  to  a  larger  value.  This  is  illustrated  for  X  *  1,5 
by  the  lower  curve  in  figure  5.  From  the  work  that  has  been  done,  it 
appears  that  this  wave  exists  for  all  values  of  X.  This  is  contrasted 
by  the  ordinary  surface  wave  which  exists  only  for  X  greater  than  two. 

It  is  anticipated  that  a  further  study  of  the  surface  wave  modes  will 
be  made. 

7.  SUMMARY 

The  problem  of  a  point  electric  dipole  moving  over  a  dispersive 
dielectric  half-space  has  been  studied.  Expressions  for  the  fields  in 
both  reference  frames  have  been  formulated  for  the  reflected  radiation 
fields  for  the  case  of  an  arbitrary,  dispersive  dielectric  and  for  the 
lateral  wave  and  surface  wave  fields  for  the  case  of  a  lossless  plasma 
dielectric.  In  the  rest  frame  of  the  source,  it  is  found  that  all  three 
waves  exhibit  the  freouency  of  the  source  and  that  In  the  rest  frame  of 
the  dielectric,  they  all  exhibit  some  form  of  Doppler  shift.  It  is 
found  that  the  criteria  for  existence  of  the  lateral  uave  is  not  modified 
by  the  relative  motion  of  the  source  and  dielectric,  but  that  the 
criteria  for  existence  of  the  surface  wave  ia  greatly  modified  by  the 
relative  motion.  In  addition,  for  veloclcies  greater  than  some  critical 
velocity,  a  new  type  of  surface  wave  cooes  into  existence. 
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The  problem  of  a  point  electric  dipole  moving  over  a  dispersive  dielectric  half¬ 
space  is  studied.  The  dipole  is  located  in  the  free  space  above  the  dielectric  and  is 
assumed  to  be  time  harmonic  in  its  rest  frame,  oriented  perpendicular  to  the  interface 
and  moving  parallel  to  it.  Previous  work  in  this  area  has  been  mostly  confined  to  the 
use  of  a  plane  wave  as  a  source  or  a  non-dispersive  dielectric. 

Solutions  for  this  problem  are  obtained  using  integral  transform  techniaues.  The 
integration  is  performed  for  the  free  space  region  by  making  an  asvmDtotic  expansion 
for  one  integration  and  then  using  the  saddle  point  approximation  for  the  remaining 
integration.  The  solution  obtained  in  this  manner  for  the  vector  potential  is  then 
used  to  generate  the  electric  and  magnetic  fields  in  both  the  rest  frame  of  the  source 
and  the  rest  frame  of  the  dielectric.  This  yields  the  reflected  radiation  fields  for 
the  case  of  an  arbitrary,  dispersive  dielectric  and  the  lateral  wave  fields  and  surface 
wave  fields  for  the  cas«  of  a  lossless  plasma. 

The  field  patterns  are  distorted  by  the  relative  motion  of  the  source  and  die¬ 
lectric.  In  the  rest  frame  of  the  source,  all  three  waves  exhibit  the  frequency  of  the 
source  and  in  the  rest  frame  of  the  dielectric,  they  all  exhibit  some  form  of  Doppler 
shift  (different  from  the  primary  wave).  The  criteria  for  existence  of  the  lateral 
wave  is  not  modified  by  the  relative  motion  of  the  source  and  dic’<.  ^tric;  but  the 
criteria  for  existence  of  the  surface  waves  is  greatly  modified  b>  the  relative  motion, 
for  velocities  greater  than  some  critical  velocity,  a  new  type  of  surface  wave  comes 
into  existence. 
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